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ABSTRACT 

We study the conversion of a neutron star to a strange star as a possible energy source for 
7-ray bursts. We use different recent models for the equation of state of neutron star matter 
and strange quark matter. We show that the total amount of energy liberated in the conversion 
is in the range (1 — 4)xl0^'^erg, one order of magnitude larger than previous estimates, and in 
agreement with the energy required to power 7-ray burst sources at cosmological distances. 
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1. Introduction 

There is now compelling evidence to suggest that a substantial fraction of all gamma-ray bursts (GRBs) 
occur at cosmological distances (red shift z ^ 1 — 3). In particular, the measured red shift z = 3.42 for 
GRB 971214 (Kulkarni et al. 1998), and z ~ 1.6 for GRB 990123 (Kulkarni et al. 1999) imphes an energy 
release of 3 x 10^'^erg and 3.4 x lO^^'erg respectively, in the 7-rays alone, assuming isotropic emission. The 
latter energy estimate could be substantially reduced if the energy emission is not isotropic, but displays a 
jet-like geometry (Dar 1998; Kulkarni et al. 1999). Models in which the burst is produced by a narrow jet 
are able to explain the complex temporal structure observed in many GRBs (Sari et al. 1997; 1999). In 
any case, a cosmological origin of GRBs leads to the conclusion of a huge energy output. Depending on 
the degree of burst beaming and on the efficiency of 7-ray production, the central engine powering these 
extraordinary events should be capable of releasing a total energy of a few 10^'^erg. 

Many cosmological models for GRBs have been proposed. Among the most popular is the merging 
of two neutron stars (or a neutron star and a black hole) in a binary system (Paczynski 1998). Recent 
results (Janka & Ruffert 1996) within this model, indicate that, even under the most favorable conditions, 
the energy provided by vD annihilation during the merger is too small by at least an order of magnitude, 
and more probably two or three orders of magnitude, to power typical GRBs at cosmological distances. An 
alternative model is the so-called "failed supernova" (Woosley 1993), or "hypernova" model (Paczynski 
1998). 

In the present work, we consider the conversion of a neutron star to a strange star (hereafter NS^SS 
conversion) as a possible central engine for GRBs. In particular, we focus on the energetics of the NS^SS 
conversion, and not on the mechanism by which 7-rays are produced. Previous estimate of the total 
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energy E'^°™ released in the NS^SS conversion (Olinto 1987; Cheng & Dai 1996) or in the conversion of a 
neutron star to hybrid star (Ma & Xie 1996) gave E^""'" ^ lO^^erg, too low to power GRBs at cosmological 
distances. These calculations did not include the various details of the neutron star and strange star 
structural properties, which go into the binding energy release considerations. Here, we present accurate 
and systematic calculations of the total energy released in the NS^SS conversion using different models for 
the equation of state (EOS) of neutron star matter (NSM) and strange quark matter (SQM). We show that 
the total amount of energy liberated in the conversion is in the range E'^""'" =(1 — 4)xl0^^erg, in agreement 
with the energy required to power 7-ray burst sources at cosmological distances. 

The existence of strange stars (made up of degenerate u, d, and s quarks in equilibrium with respect 

to the weak interactions) is allowable within uncertainties inherent in our present theoretical understanding 
of the physics of strongly interacting matter (Bodmer 1971; Witten 1984; Farhi & Jaffe 1984). Thus, 
strange stars may exist in the universe, but until now, these have remained purely speculative entities. This 
situation changed in the last few years, thanks to the large amount of new observational data collected 
by the new generation of X-ray satellites. In fact, recent studies have shown that the compact objects 
associated with the X-ray bursters GRO J1744-28 (Cheng et al. 1998), SAX J1808.4-3658 (Li et al. 1999a), 
and with the X-ray pulsar Her X-1 (Dey et al. 1998) are good strange star candidates. Recently, Li et al. 
(1999b) have shown that the observed high and low frequency Quasi Periodic Oscillations (QPOs) in the 
atoll source 4U 1728-34 (Mcndcz & van dcr Klis 1999), are more consistent with a strange star compared 
to a neutron star, if the model of Oshcrovich & Titarchuk (1999) correctly interprets the QPO phenomena. 

Originally, the idea that GRBs could be powered by the conversion of a neutron star to a strange star 
was proposed by Alcock et al. (1986) (see also Olinto (1987)), and recently reconsidered by other authors 
(Cheng & Dai 1996). A similar model has been discussed by Ma and Xie (1996) for the conversion of a 
neutron star to a so-called hybrid star (a neutron star with a strange quark matter core). 

A number of different mechanisms have been proposed for the NS^SS conversion. All of them arc 
based on the formation of a "seed" of SQM inside the neutron star. For example, (i) a seed of SQM enters 
in a NS and converts it to a SS (Olinto 1987). These seeds of SQM, according to Witten (1984), are relics of 
the primordial quark-hadron phase transition microseconds after the Big Bang, (ii) a seed of SQM forms in 
the core of a neutron star as a result of a phase transition from neutron star matter to deconfincd strange 
quark matter (NSM^SQM phase transition). This could possibly happen when a neutron star (NS) is a 
member of a binary stellar system. The NS accretes matter from the companion star. The central density 
of the NS increases and it may overcome the critical density for the NSM— >SQM phase transition. The 
NS is then converted to a SS. In the case of accretion induced conversion in a binary stellar system, the 
conversion rate has been estimated (Cheng & Dai 1996) to be in the range (3-30) x 10^^" conversions per 
day per galaxy. This rate is consistent with the observed GRBs rate. 

However, once there is a seed of SQM inside a neutron star, it is possible to calculate the rate of 
growth (Olinto 1987; Horvath & Benvenuto 1988). The SQM front absorbs neutrons, protons, and hyperons 
(if present), liberating their constituent quarks. Weak equilibrium is then re-established by the weak 
interactions. As shown by Horvath and Benvenuto (1988), the conversion of the whole star will then occur 
in a very short time (detonation mode), in the range 1 ms - 1 s, which is in agreement with the typical 
observed duration of GRBs. A detailed simulation of the conversion process is still lacking, and only rough 
estimates of the total energy liberated in the conversion have been made. 

As wc show below, the dominant contribution to E'^""'^ arises from the internal energy released in 
the conversion, i.e. in the NSM^SQM phase transition. Moreover, the gravitational mass of the star will 
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change during the conversion process, even under the assumption that the total number of baryons in the 
star is conserved. 

The total energy released in the NS^SS conversion is given by the difference between the total binding 
energy of the strange star BE(SS) and the total binding energy of the neutron star BE(NS) 

^conv ^ BE{SS) - BE{NS) (1) 

In the present work we assume that the baryonic mass Mb of the compact object is conserved in the 
conversion process, i.e Mb{SS) — Mb{NS) = Mb- Then E'^°™ is given in terms of the difference between 
the gravitational mass of the NS and SS: = [Mg{NS) - Mg{SS)]c^. 

In general the total binding energy for a compact object can be written BE = BE] + BEq = 
{Mb — Mp)(? + {Mp — Mg)(?-, where BEj and BEq denote the internal and gravitational binding energies 
respectively, and Mp is the proper mass of the compact object defined as 
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P{r) , (2) 



p(r) being the total mass-energy density and m{r) the gravitational mass enclosed within a spherical 
volume of radius r. The proper mass is equal to the sum of the mass elements on the whole volume of the 
star; it includes the contributions of rest mass and internal energy (kinetic and interactions (other than 
gravitation)) of the constituents of the star. 

The total conversion energy can then be written as the sum of two contributions 

■^conv j^conv _j_ j^conv 

related to the internal and gravitational energy changes in the conversion. These two contributions can be 
written as: 

^conv ^ BEi{SS) - BEi{NS) = [MpiNS) - Mp{SS)]c^, (4) 

E'^?''" = BEg{SS) - BEg{NS) = [Mp{SS) - Mg{SS) - Mp{NS) + Mg{NS)]c^, (5) 

and these can be evaluated solving the structural equations for non-rotating compact objects (Oppenheimer 
& Volkoff 1939). To highlight the dependence of E'^""''" upon the present uncertainties in the microphysics, 
we employed different models for the EOS of both NSM and SQM. 

Recently, a microscopic EOS of dense stellar matter has been calculated by Baldo et al. (1997), and 
used to compute the structure of static (Baldo et al. 1997) as well as rapidly rotating neutron stars (Datta 
et al. 1998). In this model for the EOS, the neutron star core is composed of asymmetric nuclear matter 
in equilibrium, with respect to the weak interactions, with electrons and muons (/3-stable matter). In 
particular, we consider their EOS based on the Argonne vlA nucleon-nucleon interaction implemented by 
nuclear three-body forces (hereafter BBBl EOS). 

At the high densities expected in the core of a neutron star, additional baryonic states besides the 
neutron and the proton may be present, including the hyperons A, S, S, fl, and the isospin 3/2 nucleon 
resonance A. The equation of state of this hyperonic matter is traditionally investigated in the framework 
of Lagrangian field theory in the mean field approximation (Glcndenning 1985; Schaffncr & Mishustin 1996; 
Prakash et al. 1997). According to this model, the onset for hyperon formation in /3-stable-charged neutral 
dense matter is about 2-3 times the normal nuclear matter density (no = 0.17 fm~^). The latter result has 
been confirmed by recent microscopic calculations based on the Brueckner-Hartree-Fock theory (Baldo et 
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al. 1998). The appearance of hyperons, in general, gives a softening of the EOS with respect to the pure 
nucleonic case. In the present work we considered one of the EOS for hyperonic matter given in Prakash et 
al. (1997). 

For strange quark matter, we consider a simple EOS (Farhi and Jaffe 1984) based on the MIT bag 
model for hadrons. We begin with the case of massless non-interacting {ac = 0) quarks and with a bag 
constant B = 60 MeV/fm^: we denote the corresponding EOS as i360o. Next we consider a finite value 
for the mass of the strange quark within the same MIT bag model EOS. We take iris = 200 MeV (and 
m„ = md = 0, B = 60 MeV/fm"^, and ac = 0; hereafter EOS i3602oo)- To investigate the effect of the bag 
constant on the energy released in the NS^SS conversion, we take (almost) the largest possible value of B 
for which SQM is still the ground state of strongly interacting matter, according to the so-called strange 
matter hypothesis (Witten 1984). For massless non-interacting quarks this gives B = 90 MeV/fm'^; we 
denote the corresponding EOS as B90o. 

Recently, Dey et al. (1998) derived an EOS for SQM using a different quark model with respect 
to the MIT bag model. The EOS by Dey et al. has asymptotic freedom built in, shows confinement at 
zero baryon density, deconfincmcnt at high density, and, for an appropriate choice of the EOS parameters 
entering in the model, gives absolutely stable SQM according to the strange matter hypothesis. In this 
model, the quark interaction is described by a screened inter-quark vector potential originating from 
gluon exchange, and by a density-dependent scalar potential which restores the chiral symmetry at high 
density. The density-dependent scalar potential arises from the density dependence of the in-medium 
effective quark masses Mg, which, in the model by Dey et al. (1998), are taken to depend upon the baryon 
number density ub according to Mg = nig + 310MeV x sech{y^) , where no is the normal nuclear matter 
density, q{= u,d,s) is the flavor index, and is a parameter. The effective quark mass Mg(ns) goes from 
its constituent masses at zero density, to its current mass m,, as goes to infinity. Here we consider 
two different parameterizations of the EOS by Dey et al., which correspond to a different choice for the 
parameter v. The equation of state SSI (SS2) corresponds to u = 0.333 {v = 0.286). 



2. Results and conclusions 

To begin with, we fix as a "standard" EOS for neutron star matter the BBBl EOS (Baldo et al. 
1997), to explore how the energy budget in the NS— >SS conversion depends on the details of the EOS for 
strange quark matter. First we consider the B60q equation of state. The NS— >SS conversion based on 
this couple of EOSs, will be referred to as the BBBl^ B6Qo conversion model. Similar notation will be 
employed according to the EOS of NSM and SQM. The total conversion energy, together with the partial 
contributions, is shown in the upper panel of Fig. 1. As we can see, for Mb larger than ~ 1 Mq [i.e. 
values of the baryonic mass compatible with the measured neutron star gravitational masses) the total 
energy released in the NS^SS conversion is in the range (1 3) x lO'^'^crg, one order of magnitude larger than 
previous estimates (Olinto 1987; Cheng & Dai 1996; Ma & Xie 1996). Moreover, contrary to a common 
expectation, the gravitational conversion energy Eq^"" is negative for this couple of EOSs. To make a more 
quantitative analysis, we consider a neutron star with a baryonic mass Mb = 1.574M0 (see tab. 1), which 
has a gravitational mass Mq = 1.4O9M0, a radius R(NS) = 11.0 km, and a gravitational binding energy 
BEg(NS) = 4.497 xlO'"'"' erg. After conversion, the corresponding strange star has Mg = 1.254M0, R(SS) 
= 10.5 km, and BEg(SS) = 3.061 xlO^^ erg. The NS— >SS conversion is energetically possible in this case, 
thanks to the large amount of (internal) energy liberated in the NSM— >SQM phase transition. 
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Similar qualitative results for the total conversion energy are obtained for other choices of the two 
EOSs, but as we show below the magnitude of the two partial contributions are strongly dependent on the 
underlying EOS for NSM and SQM. The total conversion energy for the BBBl^ i3602oo model is plotted 
in the lower left panel of Fig. 1. Comparing with the previous case, we notice that the strange quark mass 
produces a large modification of the conversion energy, which is reduced by a factor between 2 3 with 
respect to the case = 0. The bag constant B has also a sizeable influence on the conversion energy. 
Increasing the value of B reduces E""'^^ and strongly modifies Eq™. This can be seen comparing the results 
for the BBBl^ S60o conversion model with those in the lower right panel of Fig. 1 for the BBBl^ S90o 
model. These results arc a consequence of the sizeable effects of the strange quark mass and of the bag 
constant mainly on the internal binding energy BEj{SS) for strange stars (see e.g. Bombaci 1999). In 
fact, all strange stars configurations within the B60o EOS are self-bound objects {i.e. BEj{SS) > 0). 
Strange star configurations within the i390o (B6O200) EOS are self-bound objects up to Mq ^ 0.8 Mq 
{Mq ^1.6 Mq), to compare with the corresponding maximum gravitational mass M^ax = 1-60 Mq 

{M^ax = 1.75 Mq). 

The results depicted in the two upper panels of Fig. 2 have been obtained using the EOS of Dey et 
al. (1998) for SQM, for two different choices of the parameter v which controls the rate at which chiral 
symmetry is restored to the quark masses at high density. For v = 0.286 (SS2), chiral symmetry is broken 
up to larger densities with respect to the case v = 0.333 (SSI). The parameter v has a strong influence 
on the internal binding energy of the strange star. In fact, we found that strange stars within the SS2 
(SSI) EOS are self-bound objects up to Mq ~ 0.7 Mq {Mq ~ 1.4 Mq), to compare with the maximum 
gravitational mass M^ax = 1-33 Mq {Mmax = 1-44 Mq). This effect is the main source for the differences 
in the calculated conversion energies for the two conversion models BBBl^SSl and BBB1^SS2. 

The next step in our study is to consider a different neutron star matter EOS, which allows for the 
presence of hyperons in the neutron star core. We consider one of the EOS (hereafter Hyp) for hyperonic 
matter given in Prakash et al. (1997). In the two lower panels of Fig. 2 we plot the total conversion energy, 
together with the partial contributions, for the Hyp— > B6O0 and for the Hyp^SSl conversion models. 
These results are in qualitative agreement with those reported in the previous figures. 

In table 1, we report the conversion energy, together with the partial contributions for the conversion 
of a neutron star with a gravitational mass Mg{NS) ^ 1.4Mq and for various conversion models. 

In the present work, we considered the conversion of a neutron star to a strange star as a possible 
energy source for 7-ray bursts. Our main focus was to perform an accurate calculation of the total released 
energy compatible with our current understanding of the microphysics of strong interacting matter. We 
showed that the total amoimt of energy liberated in the NS^SS conversion is in the range (1 — 4)xl0^'^erg, 
one order of magnitude larger than previous estimates, and in agreement with the energy required to power 
7-ray burst sources at cosmological distances. 

One of the authors (I.B.) thanks Prof. J.E. Horvath for valuable discussions during the Quark Matter 
Conference in Torino. 
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FIGURE CAPTIONS 



Fig. 1. The total energy liberated in the conversion of a neutron star to a strange star and the 
partial contributions from internal energy Ej°™ (curves labeled "Int" ) and from gravitational energy Eq^'^ 
(curves labeled "Gra") as a function of Mb- See text for details on the equations of state for neutron star 
matter and strange quark matter 



Fig. 2. Same as in figure 1, but for different conversion models 
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Table 1: Conversion to strange star of a neutron star with Mq ^ 1.4 M©, for different EOSs for NSM and 
SQM. Mb is the baryonic mass (which is conserved in the conversion process), Mg{NS) is the neutron star 
gravitational mass, and Mg{SS) is the gravitational mass of the corresponding strange star. All masses are 
in unit of the solar mass Mq = 1.989 x 10^^ g. Eq"'", Ef°^'", and E"^""^ are respectively the gravitational, 
internal, and total conversion energy, in unit of 10^^ erg. 
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